Isoscalar Giant Dipole Resonance and Nuclear Matter Incompressibility Coefficient 
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We present results of microscopic calculations of the strength function, S(E), and a-particle 
excitation cross sections cr(E) for the isoscalar giant dipole resonance (ISGDR). An accurate and a 
general method to eliminate the contributions of spurious state mixing is presented and used in the 
calculations. Our results provide a resolution to the long standing problem that the nuclear matter 
incompressibility coefficient, K, deduced from a(E) data for the ISGDR is significantly smaller than 
that deduced from data for the isoscalar giant monopole resonance (ISGMR). 
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Studies of compression modes of nuclei are of particu- 
lar interest since their strength distributions, S(E), are 
sensitive to the value of the nuclear matter incompress- 
ibility coefficient, K Over the last two decades, a 
significant amount of experimental work was carried out 
to identify strength distributions of the isoscalar giant 
monopole resonance (ISGMR) in nuclei fl. At present, 
Hartree-Fock (HF) based random-phase-approximation 
(RPA) calculations for the ISGMR reproduce the exper- 
imental data for effective interactions associated with in- 
compressibility K = 210 ± 20 McV (§. 

The study of the isoscalar giant dipole resonance (IS- 
GDR) is very important since this compression mode pro- 
vides an independent source of information on K. Early 
experimental attempts to identify the ISGDR in 208 Pb 
resulted with a value of E\ ~ 21 MeV for the centroid 
energy Similar result for E\ in 208 Pb was obtain in 

recent experiments 0. Very recent and more accurate 
data on the ISGDR obtained for a wide range of nuclei 
H seems to indicate that the experimental values for E\ 
are smaller than the corresponding HF-RPA results by 3 
- 5 MeV. 

It was first pointed out in Ref. |J that corresponding 
HF-RPA results for E\, obtained with interactions ad- 
justed to reproduce the ISGMR data, are higher than 
the experimental value by more than 3 MeV and thus 
this discrepancy between theory and experiment raises 
doubts concerning the unambiguous extraction of K from 



energies of compression modes (see also Ref. |10 



In this work we address this discrepancy between the- 
ory and experiment by examining the relation between 
S(E) and the excitation cross section a(E) of the IS- 
GDR, obtained by a-scattering. We emphasize that it 
is quite common in theoretical work on giant resonance 
to calculate S(E) for a certain scattering operator F 
whereas in the analysis of experimental data of <r{E) one 
carries out distorted-wave-Born-approximation (DWBA) 
calculations with a certain transition potential. Here we 
present results of accurate microscopic calculations for 
S(E) and for a(E) with the folding model (FM) DWBA 



with transition densities pt(r) obtained from HF-RPA 
calculations. We provide a simple explanation for the 
discrepancy between theory and experiment concerning 
the ISGDR and suggest further experiments. 

We also present an accurate projection method to 
eliminate the contributions of the spurious state mixing 
(SSM) to S(E) and p t (r) of the ISGDR obtained in HF- 
RPA calculations for an operator F. The method, which 
is based on the replacement of F with a properly mod- 
ified operator F v in the calculation of S(E) and pt(r), 
is quite general and applicable for any F and for any 
numerical method used in carrying out the RPA calcula- 
tion, such as configuration space RPA, coordinate space 
(continuum and discretized) RPA and with and without 
the addition of smearing. We note that this method was 
recently used in the continuum RPA calculation of S(E) 
in Ref. [OJ. Some preliminary results of our calculation 
were given in Ref. JT^ ]. 

In self-consistent HF-RPA calculation one starts by 
adopting specific effective nucleon-nucleon interaction, 
V\2, carries out the HF calculation for the ground state 
of the nucleus and then solves the RPA equation using 
the particle-hole (p-h) interaction V p h which corresponds 
to V 12 . The RPA Green's function G HQ is obtained 
from 

G = G (l + V ph G )- 1 , (1) 
where Go is the free p-h Green's function. For 



(2) 



the strength function and transition density are given by 
S(E) = J2 \{0\F\n)\ 2 5(E - E n ) = ±Im [Tr(fGf)] , (3) 



Pt(r,E) = 



1 



Z(r') 



—lmG(r r, E) 
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Note that the definition (0) is consistent with 



dr ' . (4) 



1 



S(E) = 



pt(r,E)f(r)dr 



(•5) 



In fully self-consistent HF-RPA calculations, the spu- 
rious state (associated with the center of mass motion) 
T = 0, L = 1 appears at E = and no SSM in the ISGDR 
occurs. However, although not always stated in the lit- 
erature, actual implementations of HF-RPA are not fully 
self-consistent. One usually makes the following approxi- 
mations: (i) neglecting the two-body Coulomb and spin- 
orbit interactions in V p h, (ii) approximating momentum 
parts in V p h, (hi) limiting the p-h space in a discretized 
calculation by a cut-off energy E™ x , (iv) introducing a 
smearing parameter (i.e., a Lorenzian with T/2) and (v) 
some numerical approximations. Although the effect of 
these approximations on the centroid energies of giant 
resonances is small (less than 1 MeV), the effect on the 
ISGDR is quite serious since each of these approxima- 
tions introduces a SSM in the ISGDR. 

The hrst serious attempt to correct for the effect of the 
SSM on S(E) and pt, associated with /(r) = /3(r) = 
t 3 Yim(^), was presented in Ref. However, this 

method is not accurate and leads to strong reduction in 
the ISGDR strength at low energies. In Refs. 
the effect of the SSM on S(E) was ignored and was only 
considered with regard to the energy weighted sum rule 
(EWSR) and the derivation of the collective transition 
density. In a very recent configuration space RPA cal- 
culation for the ISGDR |l(| a method for accounting for 
the SSM was implemented for /s(r). However, the re- 
sults for S(E) and pt depends on how the smearing is 
implemented (not described in the paper). 

Let us consider scattering operators, Eq. (|^), with 

f(r) = f(r)Y 1M (Sl) , A(r) = rY 1M (Q) , (6) 
and write — ImG as sum of separable terms 

R(r',r,u) = hmG(r',r,uj) = J2pa(r)Pa(r') . (7) 



Note that p a includes the energy dependence and is writ- 
ten as 



Pa{r) = Pasir) + p al (r) 



(8) 



where Paiif) is due to SSM and p a3 , associated with the 
ISGDR, fulfills the center of mass condition (for all a) 



</iP«3> = J fi(r)p a3 (r)dr = . (9) 
/.From (0) and (§) we have in an obvious notation 



R — R33 + R31 + R13 + Rn 



(10) 



with Rij = p a i(r)p a j(r and the required S(E) and 
p t can be obtained from i? 33 using (0) and (0) with /. 



Since R33 is not known we look for a projection operator 
that projects out p a i(r), 



F v = j2^( r i) = F- V F 1 



(11) 



with f v = f — vfi ■ Using (||) and ( |To| ) we have 

S V (E) = (f. n Rf v ) = (fR 33 f) + 2(fR 31 f v ) + (f v Ruf v ). 

(12) 

Note that (/r/^11/7)) is minimum for 

n = (fRufi}/(fiRuh) = E^i></i^i>/E^i> 2 

(13) 

and for the last two terms in ( |l2| ) to vanish we must have 
if Pal) = T](fipai), for all a. (14) 

The condition ( [li]) holds in case (0) has only one term, 
such as in a discretized calculation without smearing or 
in a configuration space calculation, and/or in case all 
Pai(f) are proportional to the same transition density, 
the coherent spurious states transition density |jl8f 

Pai(r) = a a p ss (r) = a a ^-Y 1M (n) , (15) 

where po is the ground state density of the nucleus. The 
value of 77 associated with p ss is then given by 



V = (fPss)/(flPs 



(16) 



To determine p t for the ISGDR we first use (|), (0), 
|, (|) and (|lj) with F v and obtain 



Pv( r ) 



J2c a [p a3 (r) + p al (r)] , 



(17) 



with c a = (f v Pa3)- To project out the spurious term from 
( |l7| ) we make use of (^) with p al = a a p ss and obtained 

Pt(r) = p v (r) - ap ss , a = (fipn)/(hp ss ) ■ (18) 

It is important to emphasize that using /, ; in (^|) with 
Pt(r) from (p^[), one obtains the required S V (E). This 
is due to the fact that the averaging process in (0) was 
carried out using F v and not F. Use of F in (|j) may 
produce erroneous results for pt{v) in (|l§| ) in case there 
are several terms in (0) . 

We now limit our discussion to the operator F 3 = 

A 

/3( r i); where we have for (|TJ) 



5 '„2 



(19) 
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We note that the values of rj obtained for p a associated 
with single p-h transitions in the ITuo region and the RPA 
results for the spurious state p t , differ from that of jl9| ) 
by less than 20% @. 

We have carried out numerical calculations for the 
S(E), p t (r) and a(E) within the FM-DWBA-HF-RPA 
theory. We used the SL1 Skyrme interaction pp| |, which 
is associated with K = 230 MeV, and carried out HF 
calculations using a spherical box of R > 25 fm. For the 
RPA calculations we used the Green's function approach 
with mesh size Ar = 0.3 fm and p-h maximum energy of 
E^ x = 150 MeV (we include particle states with princi- 
ple quantum number up to 12), since it is well-known that 
in order to extract accurate Pt(f), E™£ x should be much 
larger than the value required (E^ x ~ 50 MeV) to re- 
cover EWSR. Since in our calculation we also neglected 
the two-body coulomb and spin-orbit interactions, the 
spurious state energies differ from by a few MeV. We 
therefore renormalized the strength of the V p h by a fac- 
tor (0.99 and 0.974 for 116 Sn and 208 Pb, respectively), to 
place the spurious state at E=0.2 MeV. We have included 
a Lorenzian smearing (r/2 = 1 MeV) and corrected for 
the SSM as described above. We carried out the FM- 
DWBA calculation for <r(E) using a density dependent 
Gaussian nucleon-a interaction with parameters adjusted 
to reproduce the elastic cross section, with po and Pt from 
HF-RPA, see Ref. @. 

Using the operator / = r 2 for the ISGMR we calcu- 
lated the corresponding S(E), for E up to 60 MeV. We 
recover 100% of the corresponding EWSR and obtained 
the values of 17.09 and 14.48 MeV for the centroid en- 
ergy of the ISGMR in 116 Sn and 208 Pb, respectively. The 
corresponding recent experimental values are 16.07±0.12 



case Eq. (]?]) has one term. Thus, in configuration RPA 
calculation of pt, one may use and correct for the SSM 
contribution before the smearing process. 



and 14.17 ± 0.28 MeV, respectively |21 . 

In Figure 1, we present the results of the ESs(E) 
and ESi(E) for 116 Sn, obtained by using Eq. (§) with 
fz and fi, respectively. Note that from (||) and (|l(i| ) 
Sx(E) = (/1.R11/1) = Y.ifiPai) 2 , provides a measure to 
the contribution of the SSM to 83(E). The large con- 
tribution from the Lorenzian tail of the spurious state is 
clearly seen in the figure. Also shown in Figure 1, is the 
ratio (fzRfi)/ (fiRfi). At low energy this ratio is very 
close to 77, reflecting the fact that the transition density 
at the spurious state energy is close to that of Eq. ( O ) 
|U| . At higher energies, this ratio exhibits fluctuations 
due to the nonnegligible terms R31, Eq. (fTo|) . 

Figure 2, exhibits the strength functions for the IS- 
GDR in 208 Pb obtained from Eqs. (|), (@) and @. The 
solid line describes the result obtained using f, r Note 
that this result coincides with S V (E), which is free of 
SSM contribution. Similarly, the dashed line describes 
the erroneous result obtained using ^3 (it is also different 
from 83(E)). We find that when using /3, the excitation 
strengths obtained for certain states are sensitive to the 
value of r. The result obtained with ^3 coincides with 
that obtained with f v for r — > 0, as expected since in this 
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FIG. 1. Energy weighted strength functions for ISGDR in 
116 Sn, obtained from Eq. (|^) for the scattering operators fz 
(thick line) and /1 (thin line) with rj — §(r 2 ) = 35.6 fm 2 . 
Also shown is the ratio (fsRfi) / (fiRfi) (dashed line). 
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FIG. 2. Strength functions for the ISGDR in 208 Pb ob- 
tained from Eqs. (^]), ( |li| ) and (|^), using f$ (dashed line) and 
f v = f 3 — rjfi (solid line), with r\ — 52.1 fm 2 . 

Our results for the ISGDR, S V (E), indicate two main 
components with the low energy component containing 
close to 30% of the EWSR (for E up to 23 and 19 MeV 
for 116 Sn and 208 Pb, respectively), in agreement with the 
experimental observation |22]. Similar results were ob- 
tained for other nuclei and in other calculations ]P2] , ^9| , ^0[ 

In Figure 3 we present results of microscopic calcu- 
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lations of the excitation cross section of the ISGDR in 
116 Sn by 240 MeV a-particle, carried out within the FM- 
DWBA. The dashed lines are obtained using p co ii(r) of 
the ISGDR P,[l7|]. It is seen from the upper panel that 
the use of p co u increases the EWSR by at least 10% and 
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FIG. 3. The ISGDR in 116 Sn. The middle panel: max- 
imun double differential cross section obtained from p t (RPA). 
The lower panel: maximun cross section obtained with p co u 
(dashed line) and pt (solid line) normalised to 100% of the 
EWSR. Upper panel: The solid and dashed lines are the ra- 
tios of the middle panel curve with the solid and dashed lines 
of the lower panel, respectively. 

may shift the centroid energy by a few percent. An im- 
portant result of our calculation is that the maximum 
cross section for the ISGDR drops below the current ex- 
perimental sensitivity of 2 mb/sr/MeV for excitation en- 
ergy above 35 MeV (30 MeV for 208 Pb), which contains 
about 20% of the EWSR. This missing strength leads 
to a reduction of more than 2.5 MeV in the ISGDR en- 
ergy and thus explains the discrepancy between theory 
and experiment. More sensitive experiments and/or with 
higher a-particle energy are thus needed. 

In summary, we described and applied an accurate and 



general method to eliminate the SSM contributions from 
S(E) and pt- Our results indicate: (i) Existence of non- 
negligible ISGDR strength at low energy and (ii) Ac- 
curate determination of the relation between S(E) and 
<j(E) resolves the long standing problem of the conflict- 
ing results obtained for K, deduced from experimental 
data a(E) for the ISGDR and data for the ISGMR. 
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by the US Department of Energy under grant no. DOE- 
FG03-93ER40773. 
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